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Air stripping designs are considered to reduce the presence of volatile organic compounds in the purified water. 
Components of the wastewater streams are ranked by Henry’s Law Constant and the suitability of air stripping in the 
purification of wastewater in terms of component removal is evaluated. Distillation processes are modeled in 
tandem with air stripping to demonstrate the potential effectiveness and utility of these methods in recycling 
wastewater on the Moon. Scaling factors for distillation and air stripping columns are presented to accoimt for the 
difference in the lunar gravitation environment. Commercially available distillation and air stripping units which are 
considered suitable for Exploration Life Support are presented. The advantages to the various designs are 
summarized with respect to water purity levels, power consumption, and processing rates. 

An evaluation of reactive distillation and air stripping is presented with regards to the reduction of volatile organic 
compounds in the contaminated water and air. Among the methods presented, an architecture is presented for the 
evaluation of the simultaneous oxidation of organics in air and water. These and other designs are presented in light 
of potential improvements in power consumptions and air and water purities for architectures which include 
catalytic activity integrated into the water processor. In particular, catalytic oxidation of organics may be useful as a 
tool to remove contaminants that more traditional distillation and/or air stripping columns may not remove. A 
review of the current leading edge at the commercial level and at the research frontier in catalytically active 
materials is presented. Themes and directions from the engineering developments in catalyst design are presented 
conceptually in light of developments in the nanoscale chemistry of a variety of catalyst materials. 


Nomenclature 

°C = Degrees Celsius 

cm = Centimeter 

CO 2 = Carbon Dioxide 

EPA = Enviro nm ental Protection Agency 
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Ft 

= 

Foot 

g 

= 

gram 

HCA 

= 

Humidity Condensate Model A 

HETP 

= 

Height Equivalent Theoretical Plate 

ISS 

= 

International Space Station 

JSC 

= 

Johnson Space Center 

kW 

= 

Kilowatt 

MCL 

= 

Maximum Concentration Level 

fig/L 

= 

Micrograms per Liter 

mL 

= 

Milliliter 

mm 

= 

Millimeter 

m/s 

= 

Meters per second 

NaOH 

= 

Sodium Hydroxide 

NASA 

= 

National Aeronautics and Space Administration 

o 2 

= 

Oxygen 

Pa 

= 

Pascal 

RD 

= 

Reactive Distillation 

VOC 

= 

Volatile Organic Compound 

W 

= 

Watt 


I. Introduction 

M any designs for water purification feature disposable parts or processes that produce waste. For example, 
carbon dioxide (C02) removal with lithium hydroxide is an effective, but not readily regenerable process [1]. 
The reaction product of C02 and lithium hydroxide is lithium carbonate, which makes this process non-regenerable. 
Biocatalytic processes are also less than favorable due to their high maintenance requirements. Distillation and air 
stripping are physico-chemical processes, which are robust and well established. This report evaluates the power 
consumptions, component removal feasibility, and scaling factors for the purification of wastewater on the lunar 
surface through these processes. 


II. Air Stripping 
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Figure 1: General Schematic of an Air Stripper [2] 

Steam strippers [3] and air strippers [4] (Figure 1) are common designs for the treatment of wastewater produced by 
industries for the removal of pharmaceutical contaminants [5], oil-related contaminants [6], VOCs in contaminated 
soil [7], and urine wastewater [8], Because catalytic oxidation [9] consumes relatively little power and can be 
effective in removing contaminants in the gas and liquid phases, air stripping and steam stripping can both be 
considered. This report takes into account air stripping. Tables M and N show the components of the combined 
wastewater streams and their Henry’s law constants when available. 

Table A: Mixture of Humidity Concentrate and Urine Wastewater. 

Note: A crew of four is anticipated to produce 6.0 kilograms (kg) of urine and 7.8 kg of humidity condensate per 
day. To this amount, 33.9 grams (g) of oxone and 14.7 g of sulfuric acid is added. It is estimated that the 
components with the light blue background may be removed with air stripping. The remaining components feature 
constants in accordance with Henry’s law, which are too high for consideration. 


HCA + Urine Wastewater Model Components 
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Table B: Mixture of Humidity Concentrate and Urine Wastewater 
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Note: A crew of four is anticipated to produce 6.0 kg of urine and 7.8 kg of humidity condensate per day. To this 
amount, 33.9 g of oxone and 14.7 g of sulfuric acid is added. The strippable components are highlighted in this 
table. 


HCA + Urine Wastewater Model Components (Cont.) 

Component 
(alcohols written 
in 

organic acids 

written in pink) 

Error! 

Bookmark not 
defined. Henry’s 
Law Constants 

Ref 

Component 

Cone. 

(HCA-Blue 

+ Urine- 

Green) 

[mg/L] 

Urine+HC 

Total 

initial 

cone. 

% 

Total 

Comp. 

Component 

Target 

Cone. 

[mg/L] 

EPA MCL 
written in 
Violet 


4-acetyl morpholine 

n/a 

n/a 

n/a 

1.092 

0.612 


0.061 

Caprolactam 

n/a 

n/a 

n/a 

11.83 

6.62 


0.066 

2-butoxyethano 

n/a 

n/a 

n/a 

0.803 


0.084 

0.00042 

Glycolic acid 

n/a 

n/a 

n/a 

10.19 

5.71 

0.60 

0.0030 

N,N- 

dimethylformamide 

n/a 

n/a 

n/a 

0.608 

0.340 


0.034 

Propionic acid 

n/a 

n/a 

n/a 

3.916 

2.19 

0.23 

0.0012 

Morpholine 




0.384 

0.215 


0.022 

2-ethoxyethanol 

n/a 

n/a 

n/a 

0.504 

0.282 

0.053 

0.00026 

Lactic acid 

n/a 

n/a 

n/a 

369 + 0.32 

159 

16.7 

0.084 

2-ethyl Hexanoic 

Acid 

n/a 

n/a 

n/a 

0.37 

0.207 


0.021 

2-(2- 

ethoxyethoxy)ethanol 

n/a 

n/a 

n/a 

0.354 

0.198 

0.037 

0.00018 

Hexanoic Acid 

n/a 

n/a 

n/a 

0.582 

0.326 

0.034 

0.00017 

1-methyl-2- 

pyrrolidinone 

n/a 

n/a 

n/a 

0.339 

0.190 


0.019 

Nonanoic acid 

n/a 

n/a 

n/a 

0.335 

0.188 


0.019 

Taurine 

n/a 

n/a 

n/a 

523 

225 

23.6 

0.12 

Histidine 

n/a 

n/a 

n/a 

116.9 

50.3 

5.3 

0.026 

L-glutamic acid 

n/a 

n/a 

n/a 

412 

177 

18.6 

0.093 

Hippuric acid 

n/a 

n/a 

n/a 

171.1 

73.6 

7.7 

0.039 

a-D-glucose 

n/a 

n/a 

n/a 

793 

341 

64 

0.32 

Creatinine 

n/a 

n/a 

n/a 

1787 

768 


0.077 

4-ethyl morpholine 

n/a 

n/a 

n/a 

2.516 

1.41 


0.014 

Urea 

n/a 

n/a 

n/a 

13400 + 

2.415 

5760 


576 

Uric acid 

n/a 

n/a 

n/a 

471 

203 

21 

0.11 


1.1 Design of an Ammonia Stripper - 99.8 Percent Removal of Ammonia 


1.1.1 Flow Rates and Column Diameter 


The modeling of the distillation of urine and humidity condensate was studied at a flow rate of 20 liters per hour 
(L/hour). For this reason, the wastewater flow rate through a stripper will be studied at the same flow rate. The 
temperature of the wastewater is estimated as 40 C. The first of the stripper designs examined was an ammonia 
stripper. Stripping has previously been used for the recovery of ammonia from urine [10]. Generally, an alkaline 
pH is preferred for ammonia stripping [11]. While the pH of the pretreated urine and humidity condensates is 
particularly acidic prior to distillation, the pH could be increased in the distillate by systematic addition of lime. At 
pH 10 and at 40 C, 95 percent of ammonia is present as the gas. These conditions allow for successful removal of 
ammonia from water with efficiency. 

In summary, a column design is specified with the dimensions of 1.2 meters (m) in height and 1.7 cm in diameter. 
The gas pressure drop is 1630 pascals per meter (Pa/m). The gas flow rate is 0.508 kilograms per square meter 
(kg/m 2 ). The liquid flow rate is 20 L/day. The stripping factor is 2. The packing is 6.35-mm diameter Raschig 
rings. 
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The packing factor for 6.35-mm ceramic Raschig rings is 1600 [12]. The dimensionless Henry’s law constant for 
ammonia at 40 C, H. is 0.001436 [46], The minimum air-to-water ratio corresponds to the condition when the 
effluent gas from the stripper is m equilibrium with the incoming water. It is represented as follows: 


fa 


\ 


/ inin 


a -a 

HC 0 


Equation (1) 


The influent concentration of a mixed stream of urine and humidity condensate, as described in Appendix A, is 21 1 
mg/L. The target concentration is 0.5 mg/L. hi the equation above, the influent concentration is C 0 and the effluent 
(or target concentration is C e ). The minimum air-to-water ratio under these conditions is calculated to be 695. 

In the next part of the design study, the cross-sectional area for the packed tower is determined. In order to calculate 
this value, the actual air-to-water ratio (taking into account the stripping factor), the air pressure drop at half the 
value for flooding, the gas loading rate, and the liquid loading rate are required. 

The air-to-water ratio was calculated with the following equation: 


fa 

[q 


\ 

/ mill 


S' 


H 


ammonia 


Equation (2) 


Where S is the stripping factor (in this case, 2) and H ammoma is the dimensionless Henry’s law constant for ammonia 
at 40 C. The air-to-water ratio is calculated to be 1390. 

Raschig rings with a diameter of 6.35 mm are reported to have a pressure drop of 4 inches of water per foot at 
flooding [13a]. The pressure drop used will be half of that value (1630 Pa/nrm). See Figure 29. 
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Figure 2: Generalized Eckert Gas Pressure Drop, Gas Loading Correlation for Random-Packed Tower, and 
Flooding for Random-Packed Towers [13b, 13c] 

In order to determine the gas loading rate, the x-value for the Eckert gas pressure drop (Figure 2) is calculated. 
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Error! Bookmark not defined. X = 0.0212 


It is estimated that the corresponding y- value to this x in Figure 29 is about 0.15 for a 2-in. water pressure drop. 
With this value in hand, the gas loading rate, G m , is calculated. 


GL = 


1 ypM-pS ' 

/^t 0.1 

C fH, 


Equation (5) 
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The water loading rate, L m , is now calculated from the relation: 

,= G - 


' a Y £t? 

yQ A Pi 


Equation (6) 


L = 0.296 


kg 


m 


With the values for the liquid flow rate and the liquid loading rate, the column area is now determined as: 

Equation (8) 


A = SPl 


A = 9.1x10“ 4 w 2 

The diameter of the column is therefore 3.4 cm. 

3.1.2 Liquid Phase Mass Transfer Coefficients of the Ammonia Stripper 

In order to determine the length of the column, the mass transfer coefficients must be determined. 
The liquid phase mass transfer coefficient is calculated from the wetted surface area of the column: 
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Equation (9) 
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Equation (11) 


Equation (11) 
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The wetted surface area for this column design 


Liquid-phase mass transfer coefficient: 
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,2/3 
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Equation (12) 


Di can be calculated through the Hayduk-Laurie correlation: 


A = 


13.26x10-’ 

{ pJ 1, ( v J 5S " 


Equation (13) 


V b is the molar volume of ammonia at the boiling point, 
cubic centimeters (cm 3 )/mole. 

The viscosity of water at 40 °C, p b is 0.563 cP. 

D, =3.11x10-’ — 


This is determined through the LeBas method to be 26.7 


5 

The liquid phase mass transfer coefficient, k b is calculated as 8.4 x 10' 5 meters per second (m/s). 

3.1.3 Gas Phase Mass Transfer Coefficients for the Ammonia Stripper 

In order to calculate the gas phase mass transfer coefficients, the gas phase diffusion coefficient for ammonia must 
be determined at 40 °C. The Wilke-Lee modification of the Hirsch-Felder-Bird-Spotz correlation is used for this: 
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Equation (15) 


Equation (16) 


Equation (17) 


9 

American Institute of Aeronautics and Astronautics 



kT 


T 


Equation (18) 


313 


' air-ammonia ^ air -ammonia 151 


= 2.07 


ee = logj 


kT 


V ** air -ammonia J 


= 0.316 


£ = (-0. 14329 - 0.48343(ee) + 0.1939(ee) 2 + 0.13612ee 3 -0.20578ee 4 + 0.083899ee 5 

Equation (20) 

% = -0.274 


/ 


kT 




= 10* =0.532 


V * air -ammonia J 


(( 


1.084-0.249 


D s = 


1 1 

1 

J 

' M v 


(T 15 \ 


1 1 


ammonia air J 


M a M 


Prir ) 2 f 

I \ ammonia-air / J 


' kT A 


V G ammonia -air J 


A =2.41 


cm 


at 0. 1 bar 


3.1.4 Gas Phase Mass Transfer Coefficient 
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3.1.5 Overall Mass Transfer Rate Constant for the Ammonia Stripper 
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3.1.6 Length of the Ammonia Stripper Column 
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It has been demonstrated that removal of ammonia may be achieved at the target levels described in the beginning of 
this report for a column, which is about 1 ft in height and 0.7 in. in diameter. This design is compatible with and 
requires basic wastewater. One architecture, which may be worth considering, is to have an ammonia stripper 
upstream from a distillation column. In this way, ammonia can be removed along with other VOCs initially and 
then organic acids (such as acetic acid) may be removed with the brine water at the distillation column. 


3.2 Design of a Methanol Stripper - 99.9 Percent Removal of Methanol 

The Henry’s law constant for methanol in 0.97 mole per kilogram (mol/kg) [58 grams per liter (g/L)] of sodium 
chloride (NaCl) has been determined experimentally to be 0.29 ± 0.04. The amount of salt in the mixed humidity 
condensate is about 1 1 g/L. hi the absence of a closer Henry’s coefficient for the salinity of the solution, this value 
is used as an approximation for the actual conditions. 

As with the ammonia stripper previously described, the liquid flow rate is defined as 20 L/hour. The temperature of 
the wastewater is estimated as 40 °C. The pH is not as great a concern for the removal of methanol so it can be 
either acidic or basic. As will be discussed later, this flexibility enables the removal of methanol or a volatile 
organic (non-organic acid) of greater or equal Henry’s law constant to be removed before or after distillation 
depending on the chosen water purification architecture. 


The initial concentration of methanol in the wastewater is 4.45 mg/L. The target concentration is 0.0042 mg/L. 


The minimum air-to-water ratio is therefore: 
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With a stripping factor, S, of 2. The air-to-water ratio is: 


( n \ 


a 


\ Q ) ^ methanol 


= 11900 
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As with the methanol stripper, Raschig rings with a diameter of 6.35 mm are chosen for the methanol stripper. The 
same gas pressure drop is used (1630 Pa/nfm). 


In order to determine the gas loadmg rate, the x-value for the Eckert gas pressure drop is calculated. 
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Error! Bookmark not defined, x = 0.00248 

It is estimated that the corresponding y-value to this x in Figure 29 would be about 0.17 for a 2-in. water pressure 
drop. With this value in hand, the gas loading rate, G m , is calculated. 
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The water loading rate, L m , is now calculated from the relation: 
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With the values for the liquid flow rate and the liquid loading rate, the column area is now detennmed as: 

A^SPl 

4 

A = 5.4x 10“ 4 w 2 

The diameter of the column is therefore 2.6 cm. 

3.2.1 Liquid Phase Mass Transfer Coefficients of the Methanol Stripper 

In order to determine the length of the column, the mass transfer coefficients must be determined. 

The liquid phase mass transfer coefficient is calculated from the wetted surface area of the column and the liquid 
diffusivity, Di. 
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Liquid-phase mass transfer coefficient 


D] can be calculated through the Hayduk-Laurie correlation: 

13.26x 10“ s 

' = ( mJ , 4 (^ r 9 


V b is the molar volume of ammonia at the boiling pomt. This is determined through the LeBas method to be 39.5 
cm 3 /mole. 
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The viscosity of water at 40 C, pi, is 0.563 cP. 


D, =2.93x10' 


k, = 0.005 1 — — 

P,D, l//,g 


The liquid phase mass transfer coefficient, k ls is calculated as 8.47 x 10' 6 m/s. 


3.2.2 Gas-Phase Mass Transfer Coefficients for the Methanol Stripper 

In order to calculate the gas-phase mass transfer coefficients, the gas-phase diffusion coefficient for ammonia must 
be determined at 40 C. The Wilke-Lee modification of the Hirsch-Felder-Bird-Spotz correlation is used for this: 


1.084-0.249 


+ — *- ( T 15 ) — 

^ methanol M a„ J V M : 


M 

methanol air 


methanol -air , 


'methanol -air 


Molecular separation at collision for methanol: 


=1.1 8 (v„J ls =1.18 0.0395 = 0.402 • nm 

1 mol ) 


Molecular separation at collision for air: 


r air = 0.3711 • nm 


Molecular separation at collision for ammonia and air: 

r — - -*=i (r —- +rJ 


methanol -air 


= 0.387 • nm 


Energy of molecular attraction for methanol: 


'methanol _ j 21*7^ 


boiling, methanol 
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Boltzman’s constant 


k = 1.3804 xl0~ 16 &' cm 
s 2 • K 

^’boiling, methanol “ 337.8 ■ K 


C 

methanol _ 

k 

Energy of molecular attraction for air: 


-^ = 78.6 
k 


0 air -methanol _ ^ methanol ^ ^ air _ /(409)f78 6) — 1 79 

k ~\k k K ' ’ 


Collision function: 

kT T 


313 


£ air-methanol S air -methanol 1 79 

k 


= 1.75 


A 


£?£? = log. 


kT 




\ G air -methanol J 


= 0.243 


| = (- 0.14329 - 0.48343(ee) + 0.1939(ee) 2 + 0. 13612ee 3 - 0.20578ee 4 + 0.083899^ 5 
| = -0.274 - 0. 1 17 + 0.01 145 + 0.001953 - 0.0007172 + 0.000071086 - 0.000002365 
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kT 
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M methanol ^ air J 
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^ methanol 



0 methanol -air ) f 
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\ a methanol -air J 


D = (( l-084-0.249V^)[5538W0^ ) =23-£ ml 
g (10 132.5)(0.387) 2 (0.419) 5 


0.01 149 \ee 
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3.2.3 Gas-Phase Mass Transfer Coefficient 


k G = 5.23 {a,D g 


a t/ j 


S ) 


V Pg D s J 


1/3 


[ a , d P Y 


jU is the air viscosity 


k G = 5.23(710 • 2.3 


0.498 


710-1.75x10" 


1.75x10 


-5 A 


1/3 


1.13-2.3x10" 


(710-0.00635) 


-2 m 

s 


k G = (5.23)(1633)(13.25)(0.407)(0.0492) = 2270 


m 


3.2.4 Overall Mass Transfer Rate Constant for the Methanol Stripper 


1 + — = 1310*5 


K r a k,a k a H 

L 1 gw 


K L a = 0.00076-5 


3.2.5 Length of the Methanol Stripper Column 


( f 

1 + 


L = 


AK L a \ S - 1 


In 


^ ^ TO J 


r „ 2.3 x 10“ 7 • m 

L = 2 x -t rln 

(5.4 x 10“ 4 )(0. 00076) 


( | 4.45 A 
0.0042 


= 7.0 • m 


3.3 An Adapted Commercial Air Stripper 

All of the companies that were contacted about air stripping of wastewater deal with much larger throughput 
volumes than are required for purification of wastewater for a four- to seven-person crew. QED Environmental 
Systems, one of the companies contacted, develops air strippers and is interested in developing an air stripper for 
these purposes. The vendor for this company has made some calculations for the power requirements for the 
removal of ethanol and ammonia from wastewater. These are shown in Figure 3. One of the features of this air 
stripper is that it is easy to clean. One problem with air stripping wastewater is the fouling of the equipment. QED 
Environmental Systems uses a sliding tray technology, so that the trays in the unit can be removed periodically for 
cleaning. While a packed column is generally preferable for smaller-scale processing, this feature indicates easy 
maintenance, which may make a tray column preferable. 
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1. Water Results 

Contaminant 

Influent 

(PPb) 

Target 

(ppb) 

4 -Tray 
Results 
(ppb) 

4-Tray 

%Removal 

6-Tray 

Results 

(ppb) 

6-Tray 

%Removal 

ammonia 

18000 

0 

6.1 

99.966 

< 1 

100.000 

ethyl alcohol 

(ethanol) 

8000 

0 

4 . 1 

99.949 

< 1 

100.000 


Air Temp: 40 C 

Flow: 0.4 1pm 

Water Temp: 40 C 

Stripper: EZ-Stacker 2.xp 
for details 

- Click Stripper Air Flow: 3.96 m3/min 


Figure 3: Ammonia Readily Removed by the Air Stripper Design by QED Environmental 


QED Environmental Systems can custom design an air stripper based on their plastic E-Z Stacker 2.XP model 
stripper. This custom design would process 5.4 liters of wastewater per hour. The blower required for operation 
consumes an average of 500 W. The temperature for operation is 40 C. The calculations were determined for 1 
atmosphere of pressure. With an influent concentration of 10 parts per million (ppm) ethanol, the system can get the 
concentration of ethanol below parts per billion (ppb) levels. A picture of the basic unit produced by QED 
Environmental Systems is provided in Figure 4. In Figure 5 is the overall design of the EZ Stacker. The dimensions 
of the custom unit are shown in Figure 6. These dimensions are for operation on the Earth surface. 
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QED EZ-Stacker Model 2.4P 


air exhaust 

4* PIPE 




Copyright QED Environmental Syttemt. Inc.. 2 CO l 


Figure 5. QED Environmental Systems EZ Stacker basic design. 

Note: The dimensions in Figure 5 do not apply to the custom design. Figure 5 presents the basic design of the unit 
which can be customized. 
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r" ^ 



58 " 



20 " 


Figure 6. Dimensions of the custom design of the QED Environmental EZ Stacker provided by Dave Fischer, Vice 
President of Technology at QED Environmental Systems, Inc. This stainless steel, four tray EZ stacker design has 
an estimated dry weight of 200 kg. 
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4.0 Distillation and Air Stripping in Tandem 



Figure 7: Outline for the Model of a Gravity-Based Purification System 

Note: The feed stream is used as cooling water for the condenser of the still. The heated feed is sent to distill. The 
distillate is then passed through an air stripper to remove light key organic compounds from the water and transfer 
them into air. The air feed would then be passed through a catalytic oxidizer to mineralize the VOCs. 

In Figure 7 is the distiller-air stripper design modeled in Aspen. The same air flow rate was chosen for this stripper 
as was calculated in the ammonia stripper discussed earlier. The 6-mm packing was selected. The suggested 
diameter for the column was suggested by Aspen to be 5.3 cm. 

The components used in this model were decided to be ethanol, acetic acid, sodium chloride, ammonia, and water. 
All of these components were added in the concentrations present in the humidity condensate wastewater. A 
simulation was run with these components in basic pH (when sodium hydroxide was added to the feed) and 
compared with a simulation run under acidic conditions (when sulfuric acid was added to the feed). The results of 
these simulation experiments are shown in Table O and Table P. 


20 

American Institute of Aeronautics and Astronautics 


Table C: Distiller Concentrations from an Air Stripper and Distiller Operating in Series from an Initially 
Alkaline Wastewater Solution 


Component Mass Fraction (Distillation) 



Feed 

Distillate 

Brine Water 

Water 

9.98 x 10' 1 

1.00 

9.89 x 10' 1 

Ethanol 

5.58 x 10' 6 

6.55 x 10* 6 

1.03 x 10‘ 7 

Acetic Acid 

1.92 x 1 0‘ 13 

3.33 x 10‘ 30 

7.76 x 10‘ 13 

h 3 o + 

1.13 x 10‘ 14 

1.54 x 10‘ 12 

1.19 x 10‘ 14 

Na + 

8.23 x 10‘ 4 

1.28 x 10‘ 30 

5.48 x 10‘ 3 

OH' 

4.23 x lO’ 4 

8.93 x 10* 6 

2.82 x 10‘ 3 

ch 3 co 2 - 

8.04 x 1 0‘ 6 

2.62 x 1 O' 29 

5.36 x 10‘ 5 

cr 

3.82 x 10' 4 

1.97 x 10- 30 

2.55 x 10‘ 3 

Ammonia 

2.11 x 10' 4 

2.39 x 10‘ 4 

1.23 x 10‘ 7 

Ammonium 

2.38 x IQ -7 

9.47 x 1 O' 6 

4.04 x 10’ 11 


Table D: Air Stripping Concentrations from an Air Stripper and Distiller Operating in Series from an 
Initially Alkaline Wastewater Solution 


Component Mass Fraction (Air Stripping) 



Distillate Feed 

Pure Water 

Dirty Air 

Water 

1.00 

1.00 

7.25 x 10‘ ; 

Ethanol 

6.55 x 10‘ 6 

1.12 x 10' 6 

3.12 x 10“ 

Acetic Acid 

3.33 x 1 O' 30 

0.00 

1.09 x 10' ; 

o 2 

0.00 

5.15 x 10' 6 

7.42 x 10‘ 

n 2 

0.00 

6.25 x 10' 7 

1.85 x 10' 

h 3 o + 

1.54 x 1 O' 12 

1.37 x 10' 10 

0.00 

Na + 

1.28 x 1 O' 30 

0.00 

0.00 

OH" 

8.93 x 10' 6 

2.58 x 1 O' 9 

0.00 

Ammonia 

2.39 x 10' 4 

2.69 x 10' 11 

1.39 x 10" 

Ammonium 

9.47 x 10' 6 

2.61 x 10' 9 

0.00 


It can be seen that at these flow rates ethanol is not readily removed from the distillate water by air stripping. The 
ammonia is removed below the target concentration. Ammonia removal is a pH-dependent process. Air stripping of 
ammonia is only effective at high pH. On the other hand, separation of ammonia from water by distillation is only 
effective at low pH. 

Since the removal of ethanol is not effective at this incoming air flow rate, the air flow rate is adjusted in Aspen as 
per the calculations for the methanol air stripper to be 13.6 times the mass flow rate of the incoming water stream. 

The next calculation was done with 0.1 -percent sulfuric acid to make the initial wastewater stream acidic (See Table 
Q). The result is that all ammonia is separated from the distillate water in the first distillation step. In the next step, 
it is observed that air stripping for the removal of ethanol is more than effective at this higher air flow rate. Acetic 
acid, on the other hand, is not effectively removed in this process. Air stripping of acetic acid is known to be 
difficult and requires huge excesses of air-to-water to achieve separation (See Table R). With this design of an air 
stripper distillation column, acetic acid can be removed only if the pH is alkaline. If the pH camiot be made alkaline 
at any stage in the process, then an additional unit will need to be added to remove the acetic (and formic) acid. 
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Table E: Distillation Concentrations from an Air Stripper and Distiller Operating in Series from an Initially 
Acidic Wastewater Solution 


Component Mass Fraction (Distillation from an imtial acidic feed) 


Distillate Brine Water Feed 


Water 

1.00 



9.87 

X 

10' 1 

9.98 

X 

10' 1 

Ethanol 

6.55 

X 

IO' 6 

1.06 

X 

10‘ 7 

5.58 

X 

10‘ 6 

Acetic Acid 

5.11 

X 

10‘ 7 

4.74 

X 

10‘ 5 

8.14 

X 

10‘ 6 

H 3 0 + 

2.31 

X 

10‘ 7 

4.03 

X 

IO’ 4 

1.08 

X 

IO' 4 

Na + 

1.28 

X 

l 0 -3° 

1.65 

X 

10‘ 3 

2.48 

X 

IO' 4 

HC1 

2.02 

X 

10-3° 

1.27 

X 

io- 9 

1.22 

X 

IO -11 

OH‘ 

5.88 

X 

10' 11 

7.49 

X 

IQ- 14 

4.39 

X 

10- 14 

CH 3 co 2 - 

7.16 

X 

10' 7 

8.93 

X 

IO’ 8 

3.57 

X 

IO' 8 

cr 

1.97 

X 

10-3° 

2.55 

X 

IO’ 3 

3.82 

X 

10‘ 4 

h 2 so 4 

5.44 

X 

10-3° 

2.57 

X 

10‘ 13 

1.21 

X 

10' 15 

hso 4 - 

5.39 

X 

10-3° 

3.12 

X 

10' 3 

2.24 

X 

IO' 4 

so 4 2 - 

5.33 

X 

10-3° 

3.44 

X 

10‘ 3 

7.57 

X 

IO' 4 

Ammonia 

6.38 

X 

IO' 29 

1.47 

X 

IQ- 10 

2.10 

X 

10‘ n 

Ammonium 

8 

.29 x 10' 

■ 27 1 

.49 x 10' 3 

2 

.23 x 10' 


Table F: Air Stripper Concentrations from an Air Stripper and Distiller Operating in Series from an Initially 
Acidic Wastewater Solution 

Note: The pH of the starting feed is dropped to 1 through the addition of sulfuric acid. The air flow rate is adjusted 
to be 13.6 times (by mass) the water flow rate in the stripper. 


Component Mass Fraction 

(Air Stripping) 


Distillate Feed 

Pure Water 

Dirty Air 

Water 

1.00 

1.00 

3.88 x IO' 2 

Ethanol 

6.55 x IO' 6 

0.00 

3.93 x 10’ 7 

Acetic Acid 

5.11 x IO’ 7 

1.73 x 1 O' 6 

9.37 x IO' 9 

h 3 o + 

2.31 x IO’ 7 

4.25 x IO’ 7 

0.00 

Na + 

1.28 x 1 O’ 30 

0.00 

0.00 

OH’ 

5.88 x 10’ 11 

8.39 x IO' 13 

0.00 

ch 3 co 2 - 

7.16 x IO’ 7 

1.32 x 1 O' 6 

0.00 

cr 

1.97 x 1 O’ 30 

0.00 

0.00 

h 2 so 4 

5.44 x 1 O' 30 

0.00 

8.27 x 1 O' 38 

hso 4 - 

5.39 x 10’ 3 ° 

0.00 

0.00 

S0 4 2 - 

5.33 x 10’ 3 ° 

0.00 

0.00 

Ammonia 

6.38 x 1 O' 29 

0.00 

5.04 x IO' 25 

Ammonium 

8.29 x 1 O' 27 

0.00 

0.00 


Beyond these basic components, there are 12 compounds for which the Henry’s law constants are not available. The 
compounds, which are not adequately removed in a smgle pass by a combination of distillation through a ten-stage 
column and air stripping are shown in Table G. 

Table G: Remaining Compounds not adequately Removed by a Combination of Air Stripping and Distillation 


Organic Acids 


Acetic Acid Removable with catalytic oxidation or distillation 

Formic Acid from alkaline solutions 


22 

American Institute of Aeronautics and Astronautics 















Glycolic Acid 


Lactic Acid 

Nitrogen containing 
Organic Compounds 

Morpholine 

Catalytic oxidation may remove this compound. 

Aprotic Organic 

Compounds 

Formaldehyde 

Catalytic oxidation may remove this compound. 

Protic Organic 

Compounds 

Ethoxyethanol 

Catalytic oxidation may remove these compounds. 

Butoxyethanol 

Ethylene Glycol 


Further investigation into air stripping may demonstrate how multiple pass processes will be useful. Multiple passes 
in air stripping could remove some organic acids. It may be useful to use a group contribution method to calculate 
the Henry’s constant for these compounds. 


III. Effects of Different Gravitational Environments on Column Size 
1.1 Tray Columns [14] 

Since both the vapor and liquid velocities change, as will be demonstrated later, there is a risk of both flooding and 
weeping when the gravitation field is changed for a given column of set dimensions. The distance between trays, 
the diameter of the perforations in the trays, heights of the weirs, and cross sectional area of the column are all 
variables, which may be taken into account when adapting the optimal conditions to the column for different 
gravitational environments. 

When considering the effect of gravity on the dynamics of a gravity-based distillation column, it is useful to think of 
the combination of the velocity of gas bubbles percolating upward through the plates and the down-coming liquid 
flowing in the opposite direction. These velocities change when the acceleration, due to gravity, changes. To begin 
with, the gas velocity, V F , is defined as Qq/A^. The formula for V F is given below: 


V r 



Equation (26) 


p L and p G are the liquid and gas phase densities and C F is the flooding coefficient. This gas velocity is directly 
related to the bubble velocity. The bubble velocity can be calculated by setting the drag forces and buoyancy forces 
equal to each other. The following equation is the equation for the drag force: 


F d = 


np L v 2 B d 2 e 


Equation (27) 


This equation is for the buoyant force: 
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F. 


( P L ~ Pa )g ml 


Equation (28) 


After setting Equations 6 and 7 equal to each other and solving for V B , the following equation is obtained: 


Error! Bookmark not defined. V B = 


4 dg 

Pl~Pg ) 

3 e 

V Pl J 


Equation (29) 


From Equation 8 it follows that 

V B oc Jg Equation (30) 

Since from Equation 5 A u =A g /V F , and since V B is proportional to V F the following equation must hold true: 


A 


n 



Equation (3 1 ) 


If the column pressure, feed, and reflux ratio are kept constant, then Q G will be the same between Earth, Mars, and 
the Moon. From here one can obtain a generalized ratio for how the area scales from one lunar or planetary body (1) 
to the next (2): 


A. 


A . 



Equation (32) 


Scaling the area of the downer is now of interest. In order to do this one must consider the liquid flow rate, QL. 
The equation for the rate of laminar flow due to gravity in a vertical pipe is as follows: 

Q l = ( a) PlS Equation (33) 

*8 Ml 

From this equation one sees the same relation as inverse proportionality between area. A, and the square root of 
acceleration due to gravity, g, as in Equation 10. See the following equation: 


Error! Bookmark not defined. A d oc — — 

vs 


Equation (34) 


Similarly to Equation 1 1, the scaling ratio is: 
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Equation (35) 


Error! Bookmark not defined. 


A 


d I 1 


A 


d 2 



Takmg Equation 9, one can determine the ratio of bubble velocities to be as follows: 




V„ 


f e, Y d l 1 


8i 

\S 2 j 


\d 2 j 


Equation (36) 


Thus, Equations 12, 15, 16 provide the scaling ratios for the plate area, downer area, and bubble velocity for a 
gravity -based column on one heavenly body to another. 

In designing a column, one can determine the size of the perforations in each plate and optimize them for the 
gravitational force experienced by the different droplets. For this reason the following relation is also of 
importance: 


Mg | _ Mg 
crD 1 <jD 


Equation (37) 


Taking into account the diameter of the perforations, we are interested in calculating how column efficiencies may 
be affected by changes in gravitational field. Equation 5 is a useful proportionality in this regard. From this 
proportionality, the ratio of the ln(E) on celestial body 1 to celestial body 2 can be determined to be: 


ln(l-£), 

Error! Bookmark not defined. 

ln(l - E) 2 


a i Y K Y *ili 


V 


V“2 J 


V "2 A K L I 2 J\ 


V B 

\ 

2 

V B 

1 z 


Equation 

(38) 


To get to the perforation diameter, one needs to consider them to have a direct influence on the bubble diameter. 
Then the bubble mass transfer coefficient can be reviewed to get a relation for how the mass transfer relates to 
changes in bubble contact time in the liquid, which is related to the bubble diameter. 


k l = 


D 


1/2 


ab 


7tt 


Equation (39) 


t is the bubble contact time in Equation 19. The contact time for the bubble in the liquid is going to be proportional 
to the size of the bubble divided by its velocity. 


K l oc 


(K\ 

v d j 


1/2 


Equation (40) 


Taking the ratio of K L between two bodies, 1 and 2, and substituting into Equation 18 gives Equation 21. 

Equation (41) 


ln(l-£|, _ 

M 

[V 


1/2 

( v ^ 

r B 2 

ln(l - E } 2 
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UJ 


v 

K B 1 J 


1.2 Randomly-Packed Columns [15] 
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Thus far, only the general case for distillation with plates has been discussed. If distillation is looked at on a small- 
scale basis, the best option is to use a fractional column with some kind of packing. 

The way Pettit (1986) determined the scaling factors in the case of random packing was to consider the liquid flow 
as laminar flow and the sum of a group of “filmlets.” Between all these filmlets is gas flowing in the opposite 
direction. From this, the volumetric flow of a given individual filmlet is: 

Qj- oc S 3 W f g cos (3 Equation (42) 

From this equation 8 is the total film thickness, W f is the filmlet width, and (3 is the filmlet angle to the direction of 
gravitational pull. Summing the filmlets according the cross-sectional area will help achieve a cross-sectional 
packmg area. To sum up the filmlets, the total width of all the filmlets must be 1A<., were 1 is the length of wetted 
packing surface per unit of cross sectional area and Ac is the cross sectional area of the column. 1 and (3 are 
independent of gravity — their values depend on the packing material. Hence, the proportionality of concern for this 
study is: 


q l = s 3 A cg 


Equation (43) 


If one makes the assumption that the filmlet thicknesses do not change when changing the acceleration due to 
gravity and if Q L is constant, the scaling ratio for change in liquid area is: 


A c i = g 2 

A 2 Si 


Equation (44) 


It is now important to scale in relation to the gas flow. In order to determine these relations, Pettit (1986) equated 
the drag force and the gravitational force as he did previously with the bubble columns. Pettit found that the 
interface for turbulent gas flow is proportional to afV 2 , where V is the gas velocity and af is the film area parallel to 
the gas flow (see Darcy friction factor). The gravitational force from the liquid is proportional to ga{8 where af8 is 
the column of liquid contained in the packed bed. Equating the drag force and gravity forces one gets this equation: 


A 

l _ 

( s 2 ) 

A 

2 

K Si j 


Equation (45) 


With this scaling factor for gas one gets that the lunar column area will increase by a factor of 2.45. If the scaling 
factor for liquid is used, the factor for scaling to the Moon will increase the column diameter by 6. 

1.3 Moon-Adapted Distillers and Air Strippers 

From Section 5.2, scaling gravitational effects according to the gas flow rate and scaling according to the liquid flow 
rate in a packed column give two different results. It may be necessary to run computational fluid dynamics 
calculations to get a clearer picture of what the requirements would be for a packed column. The data presented 
below (see Table H) describing scaled distillation columns for use on the lunar surface shows a range. The low end 
of the range is from the scaling factor of 2.45 for the gas-flow scaled column and the high end of the range is the 
scaling factor, 6, as calculated for the liquid-scaled column. The liquid scaling factor provides a more conservative 
estimate for the column design so as to avoid flooding. 
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Table H: Scaling the Distillation Column and the Air Strippers Discussed in the Text for Operation on the 
Lunar Surface 


Scaling Earth Columns for the Lunar Surface 

Distillation Column 

Diameter 

Length 

Earth 

6.4 cm 

1.0 m 

Moon 

10.0 cm- 16.0 cm 

1.0 m 

Ammonia Stripper 



Earth 

3.4 cm 

0.3 m 

Moon 

5.4 cm - 8.3 cm 

0.3 m 

Methanol Stripper 



Earth 

2.6 cm 

7.0 m 

Moon 

4.2 cm - 3.9 cm 

7.0 m 


TV. Reactive Distillation and Reactive Air Stripping 


1.1 Methods in Reactive Distillation (RD) Test and Used in Industry 

A reactive distillation or air stripping design is considered to reduce the presence of volatile organic compounds in 
the purified water. Reactive purification, in this case, integrates a reactor with a distillation column. A review of 
the literature in this field has revealed a variety of functional reactive columns in industry. The design can be 
considered for homogeneous self-catalyzed reactions such as esterification and hydrolysis, in addition to 
heterogeneous catalysis [16]. A variety of RD reaction processes have been achieved. These include ester 
hydrolysis, transesterification, metathesis, esterification, nylon synthesis, isooctane synthesis, cummulene synthesis, 
and hydration of alkenes [17] and [18]. There have been some designs which have been used for the purification of 
wastewater. In particular, acetic acid has been removed from wastewater through esterification [19]. This process 
does require the addition of a primary alcohol to the water which ads a treatment step and is less desirable for 
purification of wastewater on the lunar surface. Process intensification through the addition of a catalyst to the 
distillation column would not only likely save power and volume costs in the overall equivalent system mass (ESM) 
of the unit but may also yield better purification of water. The reason for this is that in separation limited by 
azeotrope formation under non-reactive conditions, the addition of a reactive constantly changes the concentrations 
such that the separation can proceed beyond azeotrope formation. [Integrated Chemical Processes: synthesis, 
analysis, and control] 

There have been concerns raised that when human urine is used as a fertilizer in agriculture certain pharmaceuticals 
are transferred to and pollute groundwater. In the interest of destroying these pharmaceuticals in urine ozonation 
has been sought as a pretreatment step to the use of urine as a fertilizer [20]. In particular, KMU Umweltschutz 
gmbh has developed an evaporator which includes ozonation for the purification of wastewater [21]. It has been 
demonstrated that this device will remove pharmaceutical compounds from wastewater. 
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Figure 7: Distillation/Ozonation by KMIJ Umweltschutz Gmbh 

A reactive distillation column has been designed by a group at the Slovak University of Technology (Figure 7). 
This group presented a theoretical model which it verified with a realized distillation reactor which was designed to 
remove organic chloroderivatives from wastewater [22]. 

1.2 Selection of the Catalyst for Reactive Distillation 

There are three kinds of catalytic materials, which may be used for heterogeneous catalysis: monolithic catalysts, 
membrane catalysts, and arranged catalysts. Monolithic catalysts are often honeycomb-shaped and are of a 
continuous structure. The catalytically active particles are often deposited on or inside the walls of the passages of 
this structure. Membrane catalysts are different from monolithic catalysts in that they can demonstrate selectivity in 
mass transport rates for the various compounds that may try to pass through the permeable walls of the catalyst 
substrate. The kind of catalyst, which will be considered for distillation, is the arranged catalyst. Structural 
catalysts were used in this study for reactive distillation. These catalysts are a subset for arranged catalysts. They 
are usually made of superimposed sheets of a variety of geometric arrangements, which are coated with a 
catalytically active material. 

In order to select a catalyst for reactive distillation one must be chosen from a large variety of tested compounds for 
gas-phase catalysis and wet-air oxidation catalysis since this is a two-phase system. Since in the column the 
material is vaporizing and recondensing, kinetic data available for a single-phase system may not be accurate in this 
multi -phase condition. 

Platinum and palladium on gold have also demonstrated substantial promise in catalytic oxidation of polar organics 
under mild conditions [23], Au/CeCb catalysts have been reported to achieve 100-percent formaldehyde conversion 
at 75 C [24], An 18.2 percent Mn/Al 2 0 3 catalyst with 0.1-percent palladium (Pd) has demonstrated to have 
complete combustion of methanol to carbon dioxide (C0 2 ) at 90 C [25]. 

One of the main drawbacks to wet-air oxidation (WAO) is the inability of catalysts to achieve complete 
mineralization of low molecular weight organics (such as acetic acid, propionic acid, methanol, ethanol, and 
acetaldehyde). The removal of acetic acid is usually negligible at temperatures lower than 573 Kelvin (K). A major 
benefit of WAO is that ammonia can be readily oxidized to nitrogen (N 2 ) through this process [26]. It is necessary 
for further studies to be undertaken for the catalytic oxidation with a focus on reaction under ambient pressure. The 
reactions that have been reviewed up to this time require a positive pressure of oxygen (0 2 ). This situation is not 
suitable for distillation where lower pressures are often sought to reduce the amount of heat transferred to the liquid 
in the reboiler. 
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New catalysts have been tested for the treatment of wastewater. One area of heterogeneous catalysis is dedicated to 
the destruction of halogenated organic compounds (HOCs). In this way, wastewater is detoxified by Pd on 
nanoscale supports, which catalyze the hydrodehalogenation of HOCs. Nanoscale supports previously used include 
gold, zero-valent iron, or magnetite [27] [28], 

At tins time the best catalyst for use in the reactive distillation column for wastewater purification may be the 
Au/Ce0 2 catalyst reported to achieve 100-percent formaldehyde conversion at 75 C [29], If this catalyst is 
incorporated into the distillation column, it needs to be tested for oxidation of a variety of other volatile organic 
compounds (VOCs). The durability, kinetics, and lifetime of the catalyst would also need to be determined. Along 
these lines, I contacted KMU Umweltschutz about the possibility of developing a packing material (e.g. Raschig 
rings) impregnated with a nanogold catalyst. They were interested in the concept and referred me to a colleague of 
theirs in Leipzig for further discussions. 

Thermal catalytic oxidation is preferable to photocatalytic oxidation so long as the catalyst achieved complete 
oxidation of polar organics at 100 °C or below. Photocatalytic oxidation requires an additional power source to 
operate an external light source. The design of the column with regards to the placement of the bulbs and the 
integration of the cooling water would also be a concern for the photocatalytic system. Additionally, the cooling 
water would have an additional heat load from the ultraviolet (UV) light. Since catalysts have been reported to be 
active for methanol mineralization below 100 C, the simplest solution is to use a thermally-activated catalyst. 


IV. Conclusion 

Air stripping is useful for the removal of ammonia and some of the volatile organics, which are not removed by the 
distillation process. This process requires relatively low power to remove components from water. A 
commercially-available unit was found, which would require 500 W to process 5.4 L of wastewater in one hour. 

The scaling of these purification systems for the lunar gravitational environment is taken into account. Previous 
analysis of this particular problem reveals that scaling factors can be applied to the diameter of the columns to 
achieve the same purity levels that are achieved on earth. It will also be necessary to assess the role a low gravity 
environment would play in the physics of the water/air mixture (froth) and also containment and flow. Studies in 
computational fluid dynamics may shed additional light on this process for a more precise scaling factor for 
distillation columns and air strippers operating on the lunar surface. 

The processing rate for the processes modeling in this report has been set to be 20 L/day. It may be of interest to 
determine the optimum processing time for a system architecture of an air stripper, fractional distillation column, 
catalytic oxidizer, and a C0 2 removal bed. In this way, the total processing time can be optimized with respect to 
the required blower power, heating requirements, and cooling requirements. A study comparing the benefits of the 
complete water recovery system with air stripping as compared to steam stripping may also potentially be of 
interest. 
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